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Abstract 

We study "soft leptogenesis", a new mechanism of leptogenesis which does not 
require flavour mixing among the right-handed neutrinos. Supersymmetry soft- 
breaking terms give a small mass splitting between the CP-even and CP-odd right- 
handed sneutrino states of a single generation and provide a CP-violating phase 
sufficient to generate a lepton asymmetry. The mechanism is successful if the lepton- 
violating soft bilinear coupling is unconventionally (but not unnaturally) small. The 
values of the right-handed neutrino masses predicted by soft leptogenesis can be low 
enough to evade the cosmological gravitino problem. 
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1 Introduction 



After the experimental confirmation of neutrino oscillations, leptogenesis ^ has become the 
most economical and attractive scenario to explain the cosmic baryon asymmetry. Within 
a range of neutrino mass and mixing parameters compatible with experimental data, it 
successfully reproduces the value n^/s = (0.87±0.04) x 10~^° derived from nucleosynthesis 
and CMB measurements. The see-saw mechanism |2] employed in leptogenesis requires 
the existence of right-handed neutrinos with masses close to the GUT scale. Since both 
the stability of the GUT mass hierarchy and gauge coupling unification strongly suggest 
low-energy supersymmetry, leptogenesis is more natural in a supersymmetric framework. 
Once supersymmetry is introduced, sneutrino decays offer a new channel for generating an 
asymmetry. 

In this paper we want to discuss how the sneutrino decay channel is fundamentally 
different than the neutrino channel. Supersymmetry-breaking terms remove the mass de- 
generacy between the two real sneutrino states belonging to the supermultiplet of a single 
neutrino generation |3]. They also provide a source of CP violation, and the mixing be- 
tween the two sneutrino states can generate a CP asymmetry in the decay. Although the 
scale of supersymmetry-breaking is much smaller than the right-handed neutrino mass, the 
asymmetry can be sizable because of the resonant effect jH |3] of the two nearly-degenerate 
states. Contrary to leptogenesis from neutrino decay, where at least two generations of 
right-handed neutrinos are required, a single-generation right-handed sneutrino decay is 
sufficient to generate the CP asymmetry. The soft terms, and not flavour physics, provide 
the necessary mass splitting and CP-violating phase. This new mechanism of leptogenesis, 
which we will call "soft leptogenesis" can then be an alternative or an addition to the 
traditional scenario of mixing between different flavour states. 

This paper is organized as follows. In sect. 2 we describe the one-generation see-saw 
model in presence of supersymmetry-breaking effects and compute the relevant CP asym- 
metry. In sect. 3 we rederive the asymmetry following a different field-theoretical approach, 
and comment on the effect of the initial-state coherence. The baryon-asymmetry efficiency 
factor is computed in sect. 4 by integrating the complete Boltzmann equations. Finally our 
results for the baryon asymmetry are presented and discussed in sect. 5. 

As we were completing this work, a paper has appeared [H] presenting the same idea. 

2 The CP Asymmetry 

The supersymmetric see-saw model is described by the superpotential 

W = Y,jN,L,H + ^M,,N,Nj, (1) 

where i,j = 1,2,3 are flavour indices and Ni, Li, H are the chiral superfields for the 
right-handed neutrinos, the left-handed lepton doublets and the Higgs, respectively. The 
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supersymmetry-breaking terms involving the right-handed sneutrinos Ni are 

- Csoft = mlNjN, + (^AijYijNiijH + ^Bi.MijNiN, + h.c.^ , (2) 

with standard notations. 

We will consider a single generation of because, as explained in the introduction, our 
effect survives even in this limiting case. For simplicity, we will also assume proportion- 
ality of soft trilinear terms and drop the flavour index for the coefficient A. Under these 
conditions, a CP- violating phase is still present. Indeed, with a superfleld rotation we can 
eliminate all phases from the superpotential parameters Yu and M (= Mn), and with an 
i?-rotation we can eliminate the relative phase between A and B. However, the remaining 
phase is physical. 

The right-handed neutrino has a mass M, while sneutrino and antisneutrino states 
mix in the mass matrix. Their mass eigenvectors 

N_ = ^ (^e^*/2iv - e-^*/2iVt) , (3) 

with $ = arg(i?M), have mass eigenvalues 

Ml = M^ + m^± \BM\ . (4) 

The sneutrino interaction Lagrangian in the basis of flavour (A^, N"^) and mass (A^+, A^-) 
eigenstates is, respectively, 

-Cint = N(YuHil + MY*i*H'' + AYjiH^+h.c. (5) 



—N 



Hil + {A + M)i,H +i^LN_ Hil + {A- M)iiH + h.c. (6) 



Here, for simplicity, we have set $ = choosing, from now on, a basis where A is the only 
complex parameter. 

The system of N-N"^ is completely analogous to the K^-K^ or B^-B^ system, and in 
this section we will treat it with the same formalism (see e.g. ref. 0)- Its evolution is 
determined (in the non-relativistic limit) by the Hamiltonian H = M — iV /2 where, at 
leading order in the soft terms, 

M = m{ (7) 

2M ^ / 

A* 

(8) 



^ M 

A. 1 
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Here F is the total N decay width 

r = i^^M = ^mM'. (9) 
4:71 y/2n 

With this (standard) definition, m = (yy^)ii (if)^/M sets the scale for the physical (mainly 
left-handed) neutrino masses m^, since m = Y2i under the condition ^^r? = 1. 

The eigenvectors of the Hamiltonian H are 

iVi = pN + qN^ 

= pN-qN\ (10) 

^V = $l^JEll, (11) 



Pj Mi2 - 



12 — 2 12 

We consider an initial state at t = with equal densities of and A^^ At time t, the 
state has evolved into 

N{t) = g4t)N{0) + ^g.{t)N\0) 

P 

iVt(t) = Pg4t)N{0)+g^{t)N^{0) , (12) 
Q 
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iMt^-Tt/2 



COS (AMt/2) 



•g-iAftg-rt/2 (AMt/2). 



(13) 



Here AM = M+ — M_ = and we have neglected AF with respect to AM. 
We can now compute the total integrated lepton asymmetry, defined by 
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Here / is a final state with lepton number equal to 1 and / is its conjugate. Since we 
want to exploit the enhancement due to the resonance [H E] , we will disregard any other 
subleading effects. In particular, we will neglect direct CP violation in the decay (vertex 
diagrams) and include only the effect of the N-N^ mixing (wave-function diagrams). This 
means that the decay amplitudes of the flavour sneutrino eigenstates can be immediately 
derived from the interaction Lagrangian in eq. (0), setting A = 0. We will include the 
factors cp and cb to parametrize the phase space of the fermionic (/ = Hi) and bosonic 
(/ = Hi) flnal states. Taking into account the time dependence described by eq. (|T^ . the 
CP asymmetry is given by 
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Figure 1: /S.bf, defined in eg. / TT^) . as a function of z = M/T. 



Evaluating eq. (fTT| in the limit f 12 ^ M12, we find 
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Performing the time integral 



we obtain the final expression for the CP asymmetry 

TB ImA . 



(16) 
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(19) 

cf + Cb 

It is easy to understand the origin of the different terms present in eq. (fTHjl . The factor 
AB signals the presence of supersymmetry breaking and the violation of lepton number; 
{B/M)lmA signals CP violation. The resonance effect is described by rB/{r'^ + B'^), which 
is maximal when F ~ As we move away from the resonance condition, e suffers an 
extra power suppression. 

An exact cancellation occurs between the asymmetry in the fermionic and bosonic 
channels, if Ci? = cb- Thermal effects, which break supersymmetry, remove this degeneracy. 
This happens both because of final-state Fermi blocking and Bose stimulation |H], and 
because of the effective masses acquired by particle excitations inside the plasma (for a full 
discussion of the thermal effects in leptogenesis, see ref. jHI). We find 

[1-xe- Xh)X{1, xe,XH) [1 - n^(E,)] [1 - ^^.(E^)] (20) 

(21) 

(22) 

(23) 
(24) 
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Cb = A(l, xh, xi) [1 + UBiEn)] [1 + nB{Ef)] 
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where the thermal masses for the relevant supersymmetric degrees of freedom are 

ml{T) = 2ml{T) = ^gl + ^g'y + ^-X^ (25) 

mj{T) = 2mj{T) = ^^72 + (26) 

Here g2 and gy are gauge couplings and At is the top Yukawa, renormalized at the appro- 
priate high-energy scale. The value of Asi;' as a function of z = M/T is plotted in fig. [T] 
Because of Bose stimulation, A^iT^ is positive and grows with temperature. However, for 
z < 1.2, the sum of Higgs and slepton thermal masses becomes larger than M, and the 
bosonic channel is kinematically closed. Eventually, for z < 0.8, also the fermionic channel 
becomes unaccessible. This explains the abrupt changes of Abf shown in fig. H 



3 Field-Theoretical Approach 



In this section we want to study the CP asymmetry using a different procedure. We use 
an effective field-theory approach of resummed propagators for unstable (mass eigenstate) 
particles, as described in ref. [S]. The decay amplitude /_ of the unstable external state 
iV_ into a final state / is described by a superposition of amplitudes with stable external 
states f±. Adding the contributions shown in fig. I2l we obtain 



/4iV_^ /) = /_-/+ 



in. 



where Ilab{p^) are the absorptive parts of the two-point functions for a, b 
our case, are given by 

H++ = H__ = Mr, H+_ = H_+ = -ImAT. 
The amplitude for the decay into the conjugate final state is 



Mi 



M' + zH++ 



(27) 
which, in 

(28) 
(29) 



Squaring the amplitudes and multiplying by the phase-space factors cp and cb, we 
obtain the asymmetry 
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(31) 
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The corresponding results for A^_|_ are obtained by interchanging the indices + and — . 
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Figure 2: Interfering N_ decay amplitudes for the fermionic final states. Analogous dia- 
grams exist for bosonic final states. The two-point function n+_ denoted by a blob contains 
a sum of all possible intermediate states. 



Neglecting supersymmetry-breaking in vertices, from the interaction Lagrangian in 
eq. © we obtain, up to an overall normalization, /+ = !,/_ = —« for the scalar-channel 
final state (Higgs and slepton) and /+ = !,/_=« for the fermionic channel (higgsino and 
lepton). Inserting these values in eq. (jHT|) and combining the asymmetries from and 
A^+, we obtain the final expression for the total CP asymmetry 

iTB 1mA , , , 

-Abf- (32) 



452 + r2 M 



This result agrees with eq. p8|) in the limit F ^ AM. When F ^ AM, the two 
states are not well-separated particles. Therefore, the result for the asymmetry depends on 
how the initial state is preparated. If sneutrinos (like K and B) are produced in current 
eigenstates and evolve freely {e.g. if produced in inflaton decay out of equilibrium), the 
formalism followed in sect. 2 gives the correct answer, taking into account the coherence of 
the initial state. On the other hand, if N are in a thermal bath with a thermalization time 
F~^ shorter than the oscillation time AM~^, coherence is lost and eq. (j32|l gives a more 
appropriate description. Therefore in principle we are sensitive to the details of the initial 
state. In practice, the difference is inessential since we can just recast eq. ()18|1 into eq. ()H2|1 
with a redefinition of the unknown soft parameters, A — > 2A, B 2B. In the following, 
we will use eq. (jH^ in our discussion. 



4 Solutions to the Boltzmann Equations 

The baryon asymmetry is given by 

nB__( 24 + iriH \ _J_ . . 

The first factor jJU] takes into account the reprocessing of the B—L asymmetry by sphaleron 
transitions, with the number of Higgs doublets uh equal to 2. Y^'^ = 45C(3)/(7r^(7*) is the 
sneutrino equilibrium density in units of entropy density, for temperatures much larger than 
M. For the minimal supersymmetric model with one generation of right-handed neutrinos, 
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m in eV 

Figure 3: Left: Efficiency \ri\ as a function of m for M = 10^° GeV and for two different 
initial conditions: (i) vanishing initial N abundance (solid red curve); (ii) thermal initial 
N abundance, 0) (short-dashed blue curve). Right: isocurves of \r]\ = 10*, i = 

—2, —3, —4 on the (m, M) plane for the case (i). 



the number of effective degrees of freedom is = 225. Then, we obtain 

— = -8.6 X 10~^- ri. 34 

s Abf ^ 

The efficiency factor rj describes the effects caused by: i) the sneutrino density be- 
ing smaller than the equilibrium density, ii) the wash-out from the lack of perfect out- 
of-equilibrium decay. Hi) the temperature-dependence of e through Abf- It is obtained 
by integrating the relevant Boltzmann equations. We have numerically solved the set of 
differential equations describing decay, inverse decay, and scattering processes for all su- 
persymmetric particles, including thermal masses for the particles involved p. With our 
definition of rj, the temperature-dependent part Abf has been factored out from e, see 
eq. (jnil). We have included in A^j? thermal masses and final-state statistical factors, as 
described by eqs. fEn|) - (j?T| . but we have neglected thermal corrections to the loop diagram 
generating the asymmetry (for complete expressions of the thermal corrections, see ref. 0). 

In fig. El (left) we plot the absolute value of the efficiency 77 as a function of m for fixed 
M = 10^° GeV. We consider two different initial conditions for Y^, the sneutrino density 
in units of the entropy density. In the first case, we assume that the N population is 
created by their Yukawa interactions with the thermal plasma, and set Yff{z — > 0) = 0. 
The second case corresponds to an initial N abundance equal to the thermal one, Yff{z 
0) = Y^{z — > 0). Here we are assuming that some unspecified high-energy interaction {e.g. 
GUT couplings) is responsible for bringing the sneutrinos into an equilibrium density at 
T ^ M. In fig. El (right) we present isocurves of \ri\ = 10*, i = —2, —3, —4 on the (m, M) 
plane, for the initial condition Yj^{z — > 0) = 0. This demonstrates that the efficiency is 
almost independent of M, in the range of M that is relevant for us. 
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Figure 4: Evolution of the absolute values of the abundances \Yx\ with z = M/T for 
M = 10^° GeV, m = eV (left) and m = lO'^ eV (right). Y^^\z) is denoted with 

short-dashed black line, Yfj{z) by the solid green line, while the red long-dashed line denotes 
the lepton asymmetry Yl{z) / tconst with econst = 10~^. 



The results in Fig. IHl indicate that, because of Abf, there is an extra suppression of the 
soft-leptogenesis efficiency compared to the standard leptogenesis case. Notice that this 
suppression occurs also if Yj^{z ^ 0) = Y^{z — » 0) (dashed line). The smaller m, the 
stronger the suppression, because the out-of-equilibrium decay occurs at lower T, where 
Abf is smaller, see fig. [H 

In the case Yff{z ^ 0) = (solid line) we observe a double-peak structure in \ri\. To 
understand this behaviour we plot in Fig. |3] the evolution of the abundances with z for 
M = 10^° GeV and m = 10~^ eV (left), m = 10~^ eV (right). The solid green lines denote 
Y^(z) and the red long-dashed lines denote the lepton asymmetries Yl{z)/ econst, for a fixed 
arbitrary value econst = 10~^. For reference, we also plot the equilibrium density Y^{z) 
with the short-dashed black line. 

For z < 1.2 the fermionic channel of sneutrino (inverse) decay creates an asymmetry. As 
soon as the bosonic channel is open {z > 1.2, see fig. Q), it dominates and the asymmetry 
flips sign. This is illustrated by the dip of the dashed lines at z = 1.2 in fig. 0] (both left and 
right), but this effect is inconsequential for the final asymmetry. During the A^-production 
phase, a wrong-sign asymmetry is generated compared to the right-sign asymmetry pro- 
duced in N decays. For small m (fig. |3] left) the Yukawa interactions are weak and the 
decay occurs at late time (small T) when A^f is small. Therefore the generation of the 
right-sign asymmetry cannot overcome the wrong-sign asymmetry. For larger m (fig. |3] 
right) the washout of the initial wrong-sign asymmetry is more efficient, and at late time 
an asymmetry with the right sign is created. This is observed in the right plot of fig. |3] as 
the additional sign- flip of Yl (or dip of the long-dashed curve). At an intermediate value 
of m the two effects perfectly compensate each other, and the final asymmetry vanishes, 
as shown in fig. IHl (solid line). In the case of an initial thermal A^ distribution (dashed line 
in fig. El) this cancellation never occurs, since the production phase is irrelevant. 
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5 Discussion of the Results 



We now have all the ingredients to discuss the results of the baryon asymmetry generated 
by the proposed mechanism of soft leptogenesis. The CP asymmetry is maximal when the 
parameters lie on the resonance condition, F = 2|i?|, where eq. fl32|l becomes 

e ImA , , 

(35) 



M 



From eq. (jH^ and from the results shown in fig. 13 we obtain that the presently observed 
baryon asymmetry requires^ 

M < lO^-*^ GeV. (36) 

TeV 

The resonance condition F = 2|i?| occurs when 

M^fHliXV'V « y'\o-GeV. (37) 

For typical values of B around the electroweak scale, the value of M in eq. ()37p is larger 
than what is required by eq. (jHEI), and n^/s is predicted to be too small. Soft leptogenesis 
can give a significant contribution to the baryon asymmetry only for very small values of 
B. 

Very low values of B require that the lepton-violating bilinear soft term should not 
be generated at the leading order in supersymmetry breaking, but only by some higher- 
dimensional operators. Let us consider the supersymmetry-breaking spurion superfield 
X = 6'^mMp\. Our assumption is that the leading contribution to B, coming from the 
operator J (POXMN"^ /Mpi, vanishes. In a general supergravity scenario, this is not the 
case. One can however envisage dynamical relaxation mechanisms (see e.g. ref. ^1]) which 
set i? = at leading order. Then, B is determined by the operator in the Kahler potential 
J d'^OXX^ N'^ /Mpi, which gives a value B ~ m^/M. The resonance condition F = 2\B\ in 
terms of 

Bm = VbM, (38) 



^The proportionality of the trilinear soft terms, assumed here, is certainly a questionable hypothesis 
and should not be strictly applied. However, the stability of the electroweak vacuum implies a bound on 
the size of A. Let us consider the see-saw one-generation model along a I?-flat and F/y-flat {dW/dN = 0) 
direction 



The scalar potential becomes 



2M 
Y2 



or M 



2 V 2MJ 



where x = Y (j) /M and, for simplicity, we have taken equal soft masses m for all scalar fields. Minima 
of the potential occur a,t x = [A ~ B /2 zL ^y{A — B/2)'^ — 3m^]. The request that the potential is positive 
at these minima (to avoid instabilities of the electroweak vacuum) leads to the condition \A — B/2\ < 2m. 
This shows that a departure from proportionality cannot significantly enhance the CP asymmetry, unless 
we accept to live on metastable vacua. 
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Figure 5: Regions of (m, M) plane where soft leptogenesis predicts ns/s > 0.83 x 10^^*^ for 
1mA < TeV and Bm = 100 GeV (dashed line) and 1 TeV (solid line). Soft leptogenesis 
is successful inside the contours. We have assumed a vanishing (left) or thermal (right) 
initial sneutrino density. 



IS 



M = — — — 2 X 10^ GeV. (39) 

According to our previous hypothesis, Bm is the parameter to be taken of the order of the 
electroweak scale. In this case, the value of M in eq. (j39|l is in agreement with eq. (j36|l . 

Our hypothesis of a small value of B is not technically unnatural. Indeed, 
radiative corrections to the lepton-violating bilinear term are of the form 5B ~ 
(FFl')iiAln(AVM2)/(167r2) ~ TA/M, where A is some uhraviolet cutoff scale. Thus, 
5B is much smaller than the assumed tree- level value {B ~ F). On the other hand, we 
stress that it would have been unnatural to choose a very small trilinear coefficient, since 
A receives gauge radiative corrections. 

In fig. El we quantify our results by showing the regions of parameters in the (m, M) 
plane where soft leptogenesis can predict hb/s = (0.87 ±0.04) x lO^^'' for 1mA < TeV and 
Bm between 100 GeV (dashed line) and 1 TeV (solid line). Soft leptogenesis is successful 
in the (m, M) region inside the contours. The two plots (left and right) correspond to 
vanishing and thermal initial sneutrino density, respectively. There is no overlap between 
the region of (m, M) parameters favourable for soft leptogenesis with the one suggested by 
conventional leptogenesis. 

The values of M required by soft leptogenesis (see fig. E]) are smaller than the usual 
see-saw expectation, and imply very small Yukawa couphngs, Y < 10-^(M/10^ GeVY^^. It 
should be said that soft leptogenesis is more natural in presence of a large mass hierarchy 
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of right-handed neutrinos, since one is working in the one-generation hmit. Therefore it is 
not inconsistent to predict that one generation of hes at a mass scale significant lower 
than the GUT scale. 

This result has interesting consequences for the gravitino problem. In traditional 
leptogenesis, the mass of the right-handed neutrino is bounded from below f2], M > 
2.4(0.4) X 10^ GeV for vanishing (thermal) initial neutrino densities Such values of 
M are often uncomfortably large when compared with the upper bounds on the reheat 
temperature after inflation Tjih < 10^~^°, obtained by the requirement that relic gravitinos 
do not upset the successful predictions of nucleosynthesis jH]. On the other hand, soft 
leptogenesis needs values of M in the range 10^~^ GeV, well within the limits imposed by 
the gravitino cosmological problem. 

In conclusion, we have discussed how soft leptogenesis provides an interesting in- 
terplay between lepton-number violating interactions at high energy and low-energy 
supersymmetry-breaking terms. We have found that soft leptogenesis can explain the 
observed baryon asymmetry within the range of parameters shown in fig. El This requires 
i) an unconventional (but not unnatural) choice of the lepton- violating bilinear soft param- 
eter, such that Bm in eq. (jHH|) of the order of the electroweak scale; ii) values of M in the 
range 10^~^ GeV, which is favourable to evade the gravitino problem. 
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